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ABSTRACT 



The Biot-Tolstov (B-T) exact impulse solution of diffrac- 
tion by an infinite half-plane is compared to the usual 
Helmhol tz-Kirchhof f (H-K) integral formulation and to the exact 
continuous wave (CW) solution of Macdonald. For backscatter 
the B-T and H-K solutions are found to differ significantly, 
especially near the surface of the half-plane, where the B-T 
solution gives close agreement with experiment. For forward 
scatter the two exact solutions and experimental data are in 
agreement. B-T is found to agree well with measurements of 
diffraction by a barrier perpendicular to a rigid base. By 
considering source and source image in the base separately the 
concept of "image of the source in the barrier" is found to be 
unnecessary. Use of the time domain form of the B-T solution 
in calculating the forward diffraction near a corner and be- 
hind a thin strip is shown to give results which agree well 
with measured data. Secondary diffraction effects are ob- 
served in the measurements of diffraction by a thin strip, a 
non-vertical barrier and a thick edge. 
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1 . INTRODUCTION 



Many practical problems in acoustics are complicated by 
the presence of boundaries which influence the sound field 
in some way. When the direct path between source and receiver 
is not obstructed, these boundaries represent a perturbation 
to the free field solution. If the boundaries are infinite 
in extent, the image or normal mode methods are often applied. 
Finite boundaries introduce the additional phenomenon of 
diffraction or scattering, which is often treated as an 
additional perturbation to the free-field plus reflection 
solution. When the finite boundary is such that it blocks 
the direct path the problem is somewhat different. 

In this case, the only sound reaching the receiver is 
that which has been diffracted at the edge of the finite 
boundary (assuming the boundary to be rigid) . In terms of 
airborne noise control a finite boundary that blocks the 
direct path is often called a barrier or screen. When a 
similar technique is applied to the reduction of unwanted 
sound in a shipboard sonar system, the boundary is generally 
called a baffle. The location and sice of baffles on sub- 
marines and surface ships are severely constrained by the 
size of the ships themselves. As a result, the baffle is 
often close to the radiating source and the receivers close 
to the baffle. Care must be used in selecting a theoretical 
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approach to solving the diffraction problem under these cir- 
cumstances since some of the techniques in general use involve 
simplifying assumptions which may not be valid under these 
rigorous conditions. 

Much work on the solution to the barrier problem may be 
found in the noise control literature. In this paper, two 
solutions used in noise control are examined, along with the 
little known solution by Biot and Tolstoy, and the theoreti- 
cal results are compared. In addition, experimental results 
relating to the application of theoretical solutions to 
finite barriers are presented. 
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II. THEORETICAL APPROACH 



A. GENERAL 

The usual approach to solving the problem of diffraction 
by a plate or barrier is to approximate the barrier mathe- 
matically as an infinitely thin, opaque, half-plane. This 
problem was first solved by Sommerfeld [Ref. 1] for plane 
continuous waves and then by Carslaw [Ref. 2] for spherical 
continuous waves. Macdonald's [Ref. 3] work for spherical 
continuous waves is often referenced as the rigorous solu- 
tion to diffraction by a wedge (and hence the half-plane) . 

In 1957 Biot and Tolstoy [Ref. 4] published their exact 
impulse solution using the normal coordinate formulation. 
Despite the fact that the Biot-Tolstoy impulse solution is 
in closed form and contains only simple functions, many 
authors in the noise reduction field still consider Macdonald 
as the only exact solution [Ref. 5] . 

A second general approach to the problem is Kirchhoff's 
approximate solution to the integral equation formulation 
of Huygens' Principle. This approach has been used exten- 
sively in both optics and acoustics and can be shown to give 
reasonable agreement with observed results, within certain 
limits , in spite of the fact that the basic simplifying 
assumptions are generally acknowledged to be incorrect 
[Ref. 6] . 
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There are additional approaches to this problem that will 
not be explored here. A good summary of the various solutions 
in the context of the noise control barrier can be found in 
Reference 5. 

Rigorous derivation of the three theories will not be 
attempted, however the development of each will be summarized 
in order that basic differences can be understood. 

B. BIOT -TOLSTOY FORMULATION 

M.A. Biot and I. Tolstoy published their landmark paper 
[Ref. 4] on the normal coordinate approach to wave propagation 
in infinite media in 1957. Within this paper was an applica- 
tion of the normal coordinate approach in solving the problem 
of diffraction from an infinite wedge which forms the basis 
for the theoretical results in this paper. The important 
features of their work in the present application are as 
follows : 

a. The solution is exact and in closed form. 

b. The solution can be expressed in elementary functions. 

c. The use of an impulse source provides insight into 
the contributions from various edges and a means of 
separating the diffracted from reflected and direct 
components . 

The extension of their results to other types of sources is 
straightforward . 

The normal coordinate approach is related to the normal 
mode description in acoustics which is most often encountered 
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in the solution of problems where the medium is bounded such 
as room acoustics, shallow water sound propagation, and 
speaker enclosures. In these cases the normal coordinates 
are the allowable modes of vibration (normal modes) of the 
system. The overall response of the system is the super- 
position of these discrete normal modes. Stated another way, 
an external source excites all modes at various amplitudes 
and phases and the solution then is the coherent sum of all 
these individual excitations. The essence of the normal 
coordinate approach to the solution of general diffraction 
problems is the extension of the principle of superposition 
of normal modes to infinite media. As the boundaries of the 
"enclosure" are removed toward infinity the spacing between 
normal modes becomes progressively less, with the limit 
being a continuous response spectrum with the boundaries at 
infinity. If we now introduce a source which "turns on" at 
time t = 0 (a transient source) , the solution contains only 
waves traveling outward and no radiation condition need be 
imposed. The solution summarized here is taken primarily 
from Reference 7. 

In order to formulate the problem of diffraction by a 
wedge, a cylindrical coordinate system is chosen with the 
origin at the wedge crest and the z-axis along the crest as 
shown in Figure 1. In this system the acoustic wave equation, 
in terms of displacement potential <j> , is written as 
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Figure 1, General Infinite yedge Diffraction Coordinate 
System, 
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By separation o£ variables the harmonic solutions are 



, ±iv0 u ( i , 2 ) r -> ±iy z ±iut 
cj> => e H u ’ J (Kr) e ' e 



CII-2) 



where the separation constants are related by 

,2 



2 03 

Y = - •< 



k = “ 
c 



Applying the boundary conditions for a rigid wedge, 



|f = 0 ; 0 = o, 9 = 6 



w 



restricting the amplitude to be finite at r = 0 , and choosing 
a location for our point source (r=r,0=0,z=O) that 
will result in z-axis symmetry, the solution is written as 



(j) = J (xr)cosv 0 cosy ze 

n v J n n 

n 



± io3t 



(I I - 3a) 



v 



n-rr 



n 0 



Cl I - 3b) 



w 



The solution can be separated into space - dependent and 
time dependent parts, 



4> = q 

n M n n 



(II - 4a) 



= J (xr)cosv 0cosyz 
n v n 

n 



( 1 1 - 4b) 



The q are the normal coordinates which satisfy the 
n n 



differential equation 

2 Qn 

q + 03 q = — 
n n n n u 

n 



CII-5) 
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here (3 is the generalized source function and u is pro- 
n n r 

portional to the kinetic energy in each mode. The series 
represented by (II-4) is known to be orthogonal and may be 
normalized by the condition 



U n = p 



1 0w 


f *“ r 


de 


I dz 1 



(VJj ) 2 r dr 



The medium is infinite in the r and z directions resulting 
in divergent integrals. These are handled by Tolstoy using 
a symbolic limiting procedure (discussed in detail in Re- 
ference 7, section 8-2) resulting in the following expression 



U = P 0 w 11 



(II " 6) 



n 2 w c 2 < dK dy 
It is at this point that we must choose an appropriate form 
for the generalized source, Q . 

The usual approach (and the one which we will ultimately 
use) is to choose the delta function or unit impulse source. 
It has two distinct advantages over other forms. First, it 
will result in a solution identical to that achieved via the 
Green's function approach, at least for those modes which are 
physically realizable. Second, it provides a simple building 
block with which the solution for other source forms can be 
constructed through the application of Duhamel's theorem. 

The result is analogous to analysis of electronic filters: 
the response of the system to a unit impulse input is the 
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system impulse response. The system transfer function is 
the Fourier transform of this impulse response. 

Biot and Tolstoy choose the source such that the displace- 
ment potential at a distance R is of the form 



* - - TJr Ut - R/ C) (II -7) 

where 1 (t - R /C) is the unit step function having the value 
zero for (t- R /c)<0 and one for (t- R /c)>0. This represents 
the instantaneous injection of a unit volume at the source 
coordinates and is commonly used as a mathematical approxima- 
tion to an explosion. Since the pressure and displacement 
are related by 
3 2 4 > 



P = 



8t : 



(II - 8) 



the pressure at a distance R due to the source function of 
Equation (II -8) will be 

p - 4iK3t t6Ct- R /C)] 

However, Medwin [Ref. 8] derives the solution for an impulse 
source directly from Biot-Tolstoy by simply taking the pres- 
sure to be proportional to the Biot-Tolstoy expression for 
particle velocity, since the two source functions differ by 
a first derivative. Medwin shows that this is equivalent to 
assuming a point source of strength S (volume per time) which 
flows uniformly beginning at time t = 0. The acoustic 
pressure at range R due to this source is given by 
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(II - 10 ) 



Q = - Al(-t) [V 2 </j ] 
x n 1 n o 

where X is the bulk modulus of the fluid, and if; is the space 

dependent solution for the free modes (Equation II -4). The 

Lapacian of ^ , in brackets, is evaluated at the location of 

n 

the source. Substituting (II -4) into (II -10) and utilizing 
(I I -6), the differential E quation (II -5) can be written as 

2 2 

q + w 2 q = -E— l(-t) J (icr )cosv 6 sc df<dy (II -11) 

n n n nu0 v ^v o no 

w n 

The solution to (I I - 1 1 ) represents the normal coordinates for 
the problem, 



l n 



^2^ — !£L j (Kr )cosv 0 Kdscdy t>0 

. , 2 ire,, v o J no - 

co w n 



( 11 - 12 ) 



Substituting into (II -4a) and taking the time derivative, the 
solution can be written as 



i! 

at 



2C : 

TT0 



y COSV 0 COSV 0 
L n n 



w n 



n o 



J (scr)J (icr)cosYzX 
v v o 

n n 



o o 



sin cot KdKdy (II -13) 

w 2 

Equation (11-13) represents an exact solution to the 
perturbation of a propagating wave by the presence of a 
rigid wedge, assuming Equation (II -10) as the form of the 
source. Its present form is not particularly useful, how- 
ever the application of two known integral transforms to 
Equation (11-13) reveals its simplicity and its close re- 
lationship to the physical problem in an intuitively 
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satisfying way. The application of these transforms is 
discussed next. 

Equation (11-13) contains integrals with respect to both 
the z-axis wavenumber y and the radial wavenumber k. Con- 
sidering the integration over y first, 

£t ,l Y^K 2 )‘ £ dy (II . 14) 
(y 2 +k 2 ) 2 

0 0 

and the relationship of the separation constants 
U3 2 = c 2 (y 2 + k 2 ) 

The integral on the RHS of Equation (11-14) represents a 
known integral transform, as follows [Ref. 7]: 



cosy z — — dy = — 

U ) 2 c 



,°° 



1 

J. 

c 



cosy s S - V .- Ct ^ ■ ^ ■ lc — I dy = y — J [<(c 2 t 
Cy z+ k 2 )^ 2c ° 



= 0 ct<z 



2 ) 2 ] ct>z 



(II -15) 



The existance of two distinct forms of the solution, one of 
which is zero, is a direct consequence of the choice of a 
transient source, i.e. if the source is turned on at t = 0 
there can be no effect at the receiver until a time equal to 
the source/receiver separation distance divided by the pro- 
pagation speed. In this case, z is one component of that 
separation vector. 

Using the result of (11-15) the remaining integral of 
(11-13) may be written as 
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(11-16) 



n 



J (k r ) J (kt ) J [K(c 2 t 2 -z 2 )^]Kc< 
v v o o J 

n n 



This integral has a known solution which takes on three 
different forms depending upon the relationship of the 



1 ^ 

variables r, r , and (c 2 t 2 - z 2 ) 2 . Tolstoy [Ref. 7] defines 
these three regions in terms of the propagation time of a 
pulse emitted by the source at t = 0 . The time of arrival 
(at the receiver) of the earlist possible direct pulse is 



given by 



t = — [(r-r ) 



„ 2 



1/2 



] 



i c L ^' O' 

The earliest time of arrival of the pulse which has traveled 
from source to receiver via the crest of the wedge is 



x = |[(r + rj 2 + z 2 ] 1/2 



0 C - ■ O' 

With these new variables, the three regions of the solution 
are given by 
I . t <t 



l 



I = 0 
n 



(11-17) 



II 



- i 



t <t <x I n = (^ rr 0 sinX) cosv n X 



X = arcos 



r 2 + r 2 + z 2 - c 2 t 2 



2rr 



0<X<tt 



(11-13) 



III 



x <t 
0 



_ i -v Y 

I = -(7rrr o sinh Y) sin(v n TT)e n 



Y = argcosh 



c 2 1 2 - (r 2 + r 2 + z 2 ) 
o 

2rr 



(11-19) 



18 



These three regions have a strong intuitive foundation. 

In Case I. it is obvious that a pulse transmitted by the 
source cannot be received until a finite time has elapsed, 
that time being equal to the separation divided by the pro- 
pagation speed. Case II. represents that time during which 
reflections from the face of the wedge are received. Tolstoy 
[Ref. 7] shows that when this expression for I is substituted 
into Equation (11-13), the solution takes the form of images 
of the source in the reflecting wall. Finally, Case III . is 
the diffracted wave. It is characterized by an arrival time 
which is later than the last possible image source arrival, 
as seen from Figure 1. 

The same results can also be interpreted geometrically, 
as discussed by Watson [Ref. 9]. If we rename the variables 
r, r Q , and (c 2 t 2 -z 2 )^ in (11-16) as b, c, and a, respectively, 
Watson points out that the three regions of the solution can 
be described in terms of the geometry of Figure 2 (z is taken 
to be zero for simplicity). The solution of Case II. applies 
whenever the value of side a is such that a triangle can be 
formed having sides a, b, c. It is obvious that this will 
be the case whenever a has the values 
(b-c) <a< (b+c) 

and, from the law of cosines, the angle A will be given by 

A b 2 + c 2 -a 2 rTT 9rn 

A = arccos — — (II-ZUJ 

where A can take on values between o and tt . The comparison 

with Equation (II -18) is immediate. In Case III . , a is 
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figure 2. Source, Receiver, and 7/edge Crest Geometry. 
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greater than (b+c) , the angle A becomes imaginary, and the 
expression corresponding to (11-20) must come from hyperbolic 
trigonometry : 



It can be seen from this discussion that the Biot-Tolstoy 
formulation in the time domain decomposes into direct/reflec- 
ted and diffracted components in an intuitively satisfying 
way. This is apparently due to the use of the normal coor- 
dinate formulation, which considers only modes which are 
physically realizable, and the choice of the transient form 
of the source. However, because of this distinct division 
between direct and diffracted solutions, considerable care 
must be exercised when working with a geometry where the 
arrival times of the direct and diffracted waves are close, 
such as near the boundary between the "illuminated" region 
and the geometric shadow. 

To arrive at the final solution for the diffracted wave, 
the result of Equation (11-19) is applied to the general 
solution (11-13) as follows: 

CO 

3b c , c-;~uv' 1 - 1 v a a • -vY 

— = — r — (rr smh Y) ) cosv 9 cosv 9 smv Tre n 
5t ir9 o L a no n n 



Expressing the sum in terms of exponentials, regrouping, and 
collecting in conjugate pairs, the solution can be written 
in exact, closed form as 



a u a 2 -b 2 -c 2 

A - argcosh — ^ 



(II -21) 



w 



n=0 




( 11 - 22 ) 
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ttY 



5in (J-( 5 ± 0 ±O o )) 




ib£ 2 

_ttY 

l-2e 9w cos (2— (tt± 9±9 ))+e 
'-Q ^ o J 




w 



Cl I - 23a) 



c 2 t 2 - (r 2 + r 



o 



) 



Y = argcosh 



2rr 

o 



(II -23b) 



where the term (tt± 0±9 o ) indicates that the entire expression 
within the curly brackets is the sum of four terms , each of 
which corresponds to a unique combination of (tt±9±9 0 ). 

Recalling the discussion concerning the forms of the 
impulse source, Equation (11-23) can now be used in place of 
the second partial derivative in (II -8) and the solution of 
the diffracted wave in terms of pressure can be written as 



where Y is given by Equation (II -23b) and 8 represents the 
term in curly brackets in Equation (II -23a). 

Equation (II -24) represents the diffracted pressure as a 
function of time due to a unit impulse source. To within a 
constant representing the source strength, this is analogous 
to the impulse response, h(t), of a general linear system. 
This impulse response is often Fourier transformed to obtain 
the transfer function H(w) and examined in the frequency 
domain. Since there is no known analytical transform of 



P(t) 



^-(rr sinh Y) _1 e 9w '(Sl 

AttA o 



4tt0 



(II -24) 
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of Equation (11-24), the numerical technique of computing the 
impulse response in discrete time steps and then transforming 
using the Fast Fourier Transform (FFT) is employed. To avoid 
the problem of infinite values of h(t) for the earliest 
diffracted arrival, Medwin uses an approach which computes 
the first time point by a numerical integration. All of the 
Biot-Tolstoy theoretical results presented in this paper are 
calculated in this way. 

C. HELMHOLTZ -KIRCHHOFF FORMULATION 

The Helmholtz -K irchhoff formulation is a mathematical 
statement of the heuristic description of wave propagation 
due to Huygens. Simply stated, a source at some distance 
from the surface (in this case a plate which has one straight 
edge and extends to infinity in three directions) insonifies 
all points on the surface S. Each small area element, ds , 
on the surface then acts as apoint source of spherical waves. 
At some observation point Q , the acoustic field is due to the 
incoming waves from the source directly plus the sum of the 
contributions from all of the point sources on the surface. 
Neglecting the effect of the direct incoming waves at 0, the 
field can be determined by integrating over the surface. To 
construct an exact solution in this manner one would have 
to know the amplitude and phase of each point "source" on 
the surface. This not possible, in general, and therefore 
the Kirchhoff approximation is used. The brief development 
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given herein follows that in Clay and Medwin [Ref. 10]. A 
detailed derivation can be found in most optics texts and 
in a book by Baker and Copson [Ref. 6] . 

From Green's theorem, assuming a source of time dependence 
e 1Wt , the filed at the observation point Q can be expressed 
as 



U(Q) » ^ 



3 e ’ lkR 

[U- — f- 

R 



e' lkR ( -3U^ 1 j 
) - — r— ( fif)) as 

S 



(II -25) 



where R is the distance from each element ds to the observa- 

- ikR 



tion point Q and 



is the free space Green's function. 



That the integrand is evaluated on the surface is denoted 
by the subscript s on the brackets. This known as the 
Helmholtz -K irchhoff integral. It describes the relationship 
of the field at Q to the field on S, but the problem of 
evaluating the field on S is not yet resolved. 

Equation (11-25) can be used to evaluate the diffracted 
wave field on either side of the surface. Our primary in- 
terest is the diffracted field when source and receiver are 
on opposite sides of the surface, however, when they are on 
the same side (herein referred to as backscatter) the result- 
ing expressions can be somewhat simplified. Since the purpose 
of this development is to compare the results with those of 
the normal coordinate approach, the simpler case will be 
used. For a complete treatment of this He lmholtz-K irchhoff 
integral equation approach, see Baker and Copson [Ref. 6]. 
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I f U s is defined as the incident wave field at the loca- 
tion of the surface, with the surface removed, then the 
quantity U within the integral can be approximated by 



U = MU 



(II -26) 



where M is the plane wave reflection coefficient. From 
equation (11-26) the normal derivative can be written as 



3U = _ 

an 1 3n 



(11-27) 



By Equation (11-26), the restriction has been imposed that 
the only wave field that can exist on the insonified surface 
of the barrier is that due to reflection of the incident 
wave field (i.e., there can be no scattered field on the 
surface) . The consequence of this assumption will be seen 
when the diffracted field near the surface is examined. 

The additional assumption that is required to make the 
integral in (11-25) tractable concerns the reflection co- 
efficient M. Kirchhoff's assumption considers it to be 
equal to the plane wave reflection coefficient for an infi- 
nite plane interface. Here the surface is assumed to be 
perfectly rigid and M = 1. Equation (II -25) can now be 
written as 



u(Q) - kr 



- ikR 

£ rn 

3n L s R 



) ds 



(II -28) 

When the integral in Equation (11-28) is taken over the 
surface of the half-plane, the U(Q) represents the backseat - 
tered wave field. If, on the other hand, the integral in 

Equation (11-28) is taken over the portion of the plane not 
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occupied by the surface (generally referred to as the aperture 
in optics literature), U(Q) would represent the wave field 
scattered forward into the shadow region on the side of the 
half-plane opposite the source. (Note that this approach 
would also require a change of sign due to Equation (11-27). 



D. MACDONALD'S FORMULATION 

The final theoretical development to be considered here 
is that of Macdonald [Ref. 5]. His solution, like Biot- 
Tolstoy, is an exact solution for the infinite rigid wedge, 
except that Macdonald assumes a continuous wave (cw) source. 
A derivation of his solution will not be attempted, instead, 
the results of Kawai, Fuj imoto and I tow [Ref. 11) will be 
used. The calculations in their paper are based on an ap- 
proximate form of Macdonald's solution given by Bowman and 
Senior [Ref. 12] but only the form of the exact solution is 
presented below. 

From Reference 11 the exact solution can be written in 
terms of the velocity potential as 



v = V(0 ) + V(-9 ) 
o o 



V 



ik 



f 00 (O 

Hi (T 2 + kR) 

“I YT~ dT + 

J (T 2 + 2kR) ' “ 

T R 



r 00 (i) 

Hi (T 2 + kR) 



(T 2 + 2kR) 1/2 



dT 

(11-29) 



where 

T r = sgn (tt - 9 -6 q j (kR! 



1/2 

kR) 



(II -30a) 
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1/2 



(I I -30b) 



T „ 
R 


= sgn (it - 9 — 9 ) (kR x 


- kR') 1/2 




3gn(X) - {*_ J 




R = 


(r 2 + r 2 + z 2 - 2rr 
o o 


1/2 

cos(9 +9 )) 

0 


R' = 


( r 2 + r 2 +z 2 - 2rr 
o o 


1/2 

cos (9 + 0 ) ) 

0 


Ri = 


1/2 

((r + r Q ) 2 + z 2 j 





k = co/c 

r, r Q , z, 9, 0 q as defined in Figure 1 



(II -31a) 



(I I -31b) 



(II -51c) 



and H 1 ^-® is the first order Hankel function of the first 
kind . 

The significant feature of this solution, first pointed 
out by Macdonald [Ref. 3], is that it consists of the sum of 
two velocity potentials, V(9o) and V(-9 0 ). These are iden- 
tified in the noise control literature as the source and 
the "image of the source in the barrier", respectively. This 
important interpretation of the solution will be discussed 
in a later section. 
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III. COMPARISON OF SOLUTIONS 



A. BIOT -TOLSTOY AND HELMHOLTZ -K I RCHHOFF 

Biot-Tolstov (B-T) and Helmholtz-K irchhoff (H-K) solu- 
tions will be compared by choosing a geometry in which both 
can be expressed in closed form. The backscatter case is 
chosen wherein source and receiver are in the same location. 
For this case the time domain solution for the H-K formulation 
is derived from Equation (11-28), once again following the 
development in Clay and Medwin [Ref. 10] with some minor 
changes in notation. 

Taking source and receiver to be in the same location 

(backscatter only) and assuming an impulse source, U g in 

Equation (II -28) can be written as 
- ikR 

U s - V un-n 

Making this substitution and transforming the surface 
integral into an integral along the edge [Ref. 10, pg. 323] 
Equation 11-28) can be written as follows: 



TT 
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here y, a, and R are defined as shown in Figure 3. If the 
y axis line through'} intersects the rigid half-plane, the 
integrand of equation (III -2) is represented by 



y y 

i i + 1 ] 

oo p 

(Case I ) 

If not, the integrand becomes 



[ 



r 



] 



(Case II ) 

Taking the inverse Fourier transform of (II I -2), and denoting 
the time domain solution by a small u, 

T 



u(Q) - fc 



OO - -r* 

2 - 2ikR r 

e 



■0° -TT 
? 



R : 



e lwt (j£ da 



(III -3) 



For Case I this becomes 



(t - 



«(t - ^) 

uco) = [— i £— ] + r c v 

• • 1 4ttv j l 



TT 

2 



4 vir 



4tt 



i (wt - 2kr) 



-dadf ] 



-oo-TT 

2 



(III - 4) 



Where 6 represents the Dirac delta function, a consequence 
of the choice of the Green's function of the Helmholtz equa- 
tion for the incident xvave field. 

The first bracket on the right hand side of (III -4) 
represents the reflection from one half of an infinite plane. 
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Figure 3 « Coordinate System for Backscatter from an 
Infinite Edge. 
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The second bracket contains the specular reflection from the 



other half of the plane but reduced 


by the integral expres- 


sion representing the fact that the 


plate is finite. This 


is in the form of an integral along 


the edge. Following Clay 


and Medwin [10], this is called the 


boundary wave and defined 


as follows : it 

00 2 

c 2 f f e i(ut-2kr) 

D(t) = J- J J 7 dadf 

-oo 7T r 

'I 


(I II -5) 



Case II results in a similar expression but without the 
reflection. 

In terms of Equation (III -5) the following expressions 
are written for the time domain impulse solution: 



S(t-2X) 

u(Q) - [~— C - - — D(t)] 

TT C 2 


(Case I) (III -6) 


u(Q) - D(t) 


(Case II) (III -7) 



Taking advantage of the time domain form of the solution 
and separating the reflected and diffracted (boundary wave) 
components of the wave field results in 





(-) Case I 


u(Q) = + D(t) 

TT C 2 


(+) Case II 


This can be written in terms of the 


geometry of Figure 3 as 


follows (from reference 10, section 


10.2.3) 
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u(Q) = - 



7T C ‘ 



2T 



tCt 2 + T x 2 -T 0 2 )(t : 



T0 a ) 



2,172 



t> T i 



(III -8a) 



= 0 
where 



T 



x 



V 



To 



2x 



c 

Cto 

2 

2R 

c 



t < T o 



2r 



cos 9 



To COS 0 



(III -8b) 



Finally, through trigonometric substitutions, 

u(Q) = - r 2 ! in : 9 - iy , > 

TT C 3 t (t 2 - to 2 ) (t 2 -To 2 ) ^ 



(III -9) 



The diffracted pressure at Q as a function of time is given 



bv 

p d (t, Q) = - P fi R 0 u(0) 



(III -10) 



where P. is the source pressure spectral density and R. is 
the reference distance. Taking the product P^R^ = 1 Pascal 
meter, the diffracted pressure at Q , as a function of time, 
can be written as 



P d (t, Q) = - u (Q) 

2 j- r 2 sin 20 

TTC 3 \(t 2 -t 0 2 ) (t 2 -T 0 2 ) 1/r 



(III -11) 



The corresponding expression for the Biot-Tolstoy for- 



mulation can be obtained from equation (11-22). For the 
thin plate (0 = 2ir) , 



it = £ — (rr sinh Y) ” 1 Z 
^ 2u 2 0 



Working with the sum and 
ceiver coincident), 

_nY 

„ _ 2 n0 . nir 2 

S = ZcosS^— s m— e = 

n 



nY 



n0 n0 o . nr 2 

cos cos — sin — e 



(III -12) 



noting that 0 = 0 Q (source and re 



.1 „ , e in9 - e- in V;_ nr _ 2 

Z 2 (l g s in 2 e 

n 



nY 



Multiplying terms and observing that sin — is zero for n 

even and ±1 for n odd, the sum can be evaluated as 
_ Y _3Y _ 5Y 

S = i(e 2 -e " + e 2 - ) » 

. 1 ( e ‘ u > - e' 3ui * e' 5u ' - ) + 

* 1 ( e' U2 - e" * e‘ 5u2 - ) 

4 

= j sech i + i sech ui + -g- sech u 2 ; Re(j) >0 



(HI-13) 

where 



Ul 


II 

t^|-C 
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i0 


u 2 




i0 



After considerable manipulation the sum can be written as 



y 

j 2 cosh 2 - cos 9 + cosh Y + cos 29 



s ■ T t 



(cosh Y + cos 29) cosh^- 



2 (HI-14) 

From the geometry of Figure 3 the following can be defined 

ct 0 



r - ^ 
to - — 



V = 



? 



to = 



2r ( 



2r 

c 



r = r - CT 0 

r o r - 



The various trigonometric quantities in equation (III -4) can 

now be expressed in terms of these as follows: 

sin 9 = — = ^-2- 

r o t 0 

cos 29 = 1 - 2 sin 2 9 =1-2 



to 

For r = r Q equation (II -23b) can be written as 
Y = argcosh (V) 

i r 2 r 2 r 2 

where V = (j £-*- . X) = (2 1_ - i) 

Using the half angle relations for hyperbolic functions 
leads to 



cosh | - (lev*!) 1 / 2 ) - */ tl 



Similarly , 

sinh Y = (V 2 - 1) 



2 -iH/2 



- 2^-y Ct 2 -t, 2 ) 1/2 
0 



And finally, to within a factor of 2, 



p(t) = - 



Spr ^ t 2 (cosQ + 1) 



2t 2 c 2 



t 2 (t : 



)(t 2 -x 0 2 ) 



1/2 



} = - ***o£z< 



) 

(I II -15) 



Letting P„R q be unity, as before, results in an expression 
for the diffracted pressure wave in the time domain which 
can be compared directly to the H-K result. 

Biot -Tolstoy : 

P( t ) R T= ' X V , i tt ) t>T, 



. _ 2r j- tx (cos 9 + 1) - t 0 2 

BT ’ 102 t(t 2 - t„ 2 )(t 2 - 



Helmholtz -K irchhoff : 



(III -16) 



P(t) 



HK 



zr 



7TC ' 



r sin 20 



t(t ; 



)(t : 



- x 2 ) 1/2 ' 
o 



t>x, 



(II I -17) 



Equation (III -16) and (III -17) are similar in overall 
form but differ in several important ways. When x is much 
greater than t , the time dependence of both expression is 
determined by the (t 2 -x c 2 ) ^ factor. (for t Q approaching 
x , the behavior is more complex and is treated in detail 
in Appendix A) . 

Since the impulse source has an infinite amplitude, the 
leading edge of the diffracted pressure signal (t=x ) must 
also be infinite. After this infinite leading edge, the 
time behavior can be described as follows: 

At = t - to t 2 = At 2 + 2tx 0 -To 2 

For the 'earlv diffracted signal, At<<x and 

> Of r\ 



ct 



2 - to 2 ) 1 / 2 = (2tx 0 - 2x 2 ) = (2x 0 (t -x 0 ) ) ly/ “=(2x 0 At) 



- 1/2 
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_ ^ 

The initial diffracted time signal decays as At 2 for both 

solutions. Since the amplitude is highest in the early time 

- ^ 

portion of the signal, this At 2 dominates the Fourier trans- 
form of the signal and the frequency dependance of both 

- X' 

solutions takes the form f 2 which is characteristic of many 
diffraction problems. As time increases, however, the H-K 

. - 4 

solution approaches the form t while the B-T form approaches 

- 3 

t 

Both solutions have been implemented on the IBM computer 
system at NPS. The diffracted time signals for a representa- 
tive geometry are shown in Figure 4. Care should be exer- 
cised in comparing the magnitude of the first point (t = 0) 
since the "infinite" values here were handled differently in 
the two solutions. The approach used in the B-T computation 
is the numerical integration technique discussed earlier 
while the approach for the H-K results are taken directly 
from Clay and Medwin [Ref. 10] . 

The behavior of the two solutions as a function of 0 is 
also fundamentally different. Examination of the numerators 
of Equations (III -16) and (III -17) over the region of 8 from 
0 to 180° (the problem of the thin plate being symmetrical 
with respect to 180°), reveals three apparent zero values 
for the H-K solution (9, 90°, 180°) vice two for B-T (90°, 
180°). Actually, the zero in the numerator at 90° is offset, 

in the limit, bv the term (t 2 -t 2 ) in the denominator which 

o 

goes to infinity at t = t Q . Physically, this means that 
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Kirchhoff Backseat ter 



when source and receiver are directly over the edge, perpen- 
dicular to the plate, the solution consists of reflection 
only. In order to show the detailed 8 dependence, the time 
domain solution has been computed for a variety of geometries 
and each time solution transformed to the frequency domain. 
These results are plotted versus 8 for several frequencies 
in Figures 5 and 6. The amplitudes are in dB with the 
reference being twice the range from source to edge. The 
differences between the two solutions are significant, with 
the only region of agreement being near 90° where both solu- 
tions appear to approach the value of -6dB. This is the 
expected result and represents one half of the pressure 
amplitude that would result from reflection by an infinite 
plane. Actually, the H-K formulation produces the correct 
result at 8 = 90° while the B-T solution appears to fail in 
the immediate vicinity of 90°. As 8 approaches 90°, the 
travel time of the reflection t approaches the travel time 
of the leading edge of the diffracted pulse t . 

In this region, equivalent to approaching the shadow 
boundary in the forward scattering case, the B-T solution is 
approaching the delta function form and the approximation 
used to calculate the first time point may not be valid 
(see Appendix A). The B-T diffraction solution can only be 
used for 9 approaching 90° since at 9 = 90° the separate 
reflection solution would apply (Equations 11-18 and 11-19). 
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Figure 5* Comparison of B-T and H -K Backscatter Solutions 
as a Function of 9 at r= 20 cm and a Freauency 
of 781 Hz. 
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Figure 6. Comparison of B-T and H-K Eackscatter Solutions 
as a Function of Q at r = 20 cm and a Frequency 
of 10,100 Hz. 
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As the source and receiver are moved away from 90° the 
disagreement between the solutions is dramatic. I t is in- 
tuitively correct that both solutions go to zero as 9 
approaches 180°, since the plate is infinitely thin and 
would not be "sensed” by a wave from this direction. The 
way in which the solutions approach this zero value is 
quite different. 

The most significant disagreement between the solutions 
is seen when source and receiver are close to the plate (9 
approaching 0) . The H-K solution approaches zero near the 
plate (it is symmetrical, except for sign, about 90° due to 
the sin 29 term) while B-T approaches a constant value for 
each frequency. This zero in the H-K solution appears as a 
consequence of the simplifying assumption stated in Equation 
(11-26), where the total pressure field in the plane of the 
plate is assumed to be the value of the incoming wave field 
evaluated at that location (M = 1 assumed) . This "boundary 
condition" therefore forces the diffracted pressure field 
to be zero anywhere in the plane of the plate. It will be 
shown in the experimental portion of this paper that the 
B-T solution is correct in this region near the plate. 
Although not shown in Figures 5 and 6, both solutions exhibit 
the proper sign behavior, i.e., the direct and diffracted 
pulses are opposite in sign in the region 0 to 90° but have 
the same sign for 90 to 180°. 
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The summary, the B-T and H-K solutions simplified for the 



case of coincident source and receiver, are dramatically 
different. They approach agreement, in general, only in the 
region near, but not at , 8 = 90°. 

The final consideration before leaving the topic of 
backscatter is an examination of the form of the B-T solution. 
Multiplying out the numerator and regrouping terms, Equation 
(III -16) can be written as 



r _ 2r r 
P d 1 BT ' ttc 2 *- 



cos 0 



CU-t/MU-T , 2 ) 171 t 2 (t 2 -r „ 2 ) 1/2 



) t>T 0 
(1 1 1 - 18) 



When written in this way, the solution is seen to have an 
interesting characteristic. The solution is the sum of two 
terms, one which takes into account the overall geometry of 
source/receiver and plate (the first term in III -18) and a 
term which depends on the distance from the edge of the plate 
to the source/receiver . This seems to be a characteristic 
of exact solutions of the problem of diffraction by a half- 
plane and is directly related to the "image in the barrier" 
first mentioned by Macdonald [Ref. 3] . This "image" approach 
will be discussed in some detail in the next section. 



B. BIOT -TOLSTOY AND MACDONALD 

In comparing the B-T and Macdonald half-plane diffraction 
solutions, an approach similar to that of the previous sec- 
tion will be used. However, the solutions will not be limited 
to backscatter but will be kept general and the emphasis will 
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be on "forward scatter", i.e., diffraction of sound into the 
geometric shadow region (source and receiver on opposite sides 
of the half-plane). The purpose of the comparison is twofold: 
to show that computations based on the two "exact" solutions, 
one impulse and one CW do, in fact, produce the same results, 
and to gain some insight into the use of the "image in the 
barrier". As discussed in section II, the Macdonald's 
solution form and computations are from Kawai's excellent 
paper [Ref. 11} . 

The comparison between B-T and Macdonald computational 
results requires some manipulation of geometric quantities. 
Kawai's results are converted to the notation of this paper 
for consistency, however, this results in only 2-3 points 
per frequency spectrum. To alleviate this problem somexvhat, 
the geometry is selected in such a way that several different 
B-T runs result in basically the same values. Figure 7 
shows this geometry and the comparison between B-T and 
Macdonald results. The reference pressure is the free-field 
pressure that would be measured at a distance of r + r Q from 
the source with the plate removed. Data referenced in this 
manner are labeled "dB re free field at 2r". This combination 
of geometry and reference results in B-T theoretical values 
which vary only one dB over the range 0 = 330° to 0 = 360°. 
Kawai's computations are then converted from his normalized 
presentation [Ref. 11, Figure 3] to the format of Figure 6 
in 5° increments of 0 between 330° and 360° inclusive and 
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Figure 7« Comparison of B-T and Macdonald's Forward 
Diffraction solutions for Various G's at 
r = 25 cm as a Function of Frequency. 
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plotted in Figure 7 as discrete points. The agreement between 
the two solutions is seen to be very good. The scatter in the 
Kawai points is almost certainly attributable to the process 
of reading these values from the graphs of Ref. 11. Additional 
comparisons were made over a wide range of different geometries 
using the B-T theoretical results of Bremhorst [Ref. 12] with 
agreement similar to that shown in Figure 7 observed in all 
cases . 

It should be noted, however, that the comparisons were 
not made for the general case of high frequency, large range, 
near the shadow boundary (R^ large and N small in the presen- 
tation of Kawai ) . 

Kawai [Ref. 11] actually presents two sets of computations 
of Macdonald's theory; one for the "real" source only (refer- 
red to as ATTg ) and one for the "real" plus "image" source 
o 

(ATT) . This terminology appears to be traceable directly to 
Macdonald [Ref. 3], In his solution he separates a complex 
integral [Ref. 3, p. 422] into two integrals, performs the 
integrations separately, and recombines the two into an over- 
all solution for velocity potential [Ref. 3, p. 425]. Bowman 
and Senior's approximate form of this overall solution [Ref. 

1]] is given as equations (11-29), (11-30), and (11-31). 

The comparison of B-T and Macdonald in Figure 7, is based on 
Macdonald's complete solution, ATT. Indeed, if the complete 
solution of the half-plane problem is the sum of these two 
terms, it is not clear why one would want to consider them 
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separately. Yet this is precisely what is done by many authors 
in the noise reduction field [of. Ref. 5}. 

In order to examine the "real" and "image" terms in Mac- 
donald's solution more closely, equations (11-30) and (11-31) 
are repeated here with a slight change in the way 9 and 0 o 



are expressed; 




T r = sgn (u - (0 - 9 o ))(kR 1 - kR) 1/2 


(II -30a) 


T r - = sgn (tt — (0 + 0 Q ) ) (kRi - kR') 1 ^ 2 


(II -30b) 


R - (r 2 + r 0 2 + z z - 2rr cos(0-0 ) ) 1// ~ 

0 0 


(I I - 31a) 


R" = (r 2 + r 0 2 + z 2 - 2rr cos(9 + 0 ))^ 2 

0^0 


(II -31b) 


The primed values represent the "image" source 


terms. The 



real source terms are a function of the difference between 

9 and 9 while the "image" terms are a function of their sum. 
o 

The significance of this is shown graphically in Figure 8. 
Figures 8(a) and 8(b) show two different geometries, both of 
which have the same value of the quantity 9 — 0 q . Since 
(9 — 9 c ) is only form in which the quantities 9 and 0 q enter 
real source portion of the solution, it is easily seen that 
the real source solution is independent of the orientation 
of the plate. It depends only on the locations of source 
and receiver relative to the edge. However, the geometries 
of Figures 8(a) and 8(b) are different in terms of (9 + 0 ). 

It is this term which retains the information concerning the 
orientation of the plate (half-plane) relative to source and 
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Figure 8. Geometry of the Image of the Source in the 
Barrier. 
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receiver. Since changing the sign of 0 Q is the same, geome- 
trically, as moving the source to the opposite side of the 
barrier, one can easily see why the second term of the solu- 
tion is associated with the "image of the source in the 
barrier". This interpretation is shown in Figures 8(c) and 
8(d). 

This interpretation is not unique to Macdonald's solution. 
Recalling equation (III -13), is was shown that the B-T 
solution, for the case of coincident source and receiver, 
could be written as the sum of two terms, one which depended 
only on the orientation of the edge and one which also con- 
sidered the location of the plate. Indeed, the general form 
of the B-T solution given by equation (11-24) could be 
written as the sum of a term that is a function of 9 — 0 

o 

and another that is a function of 9 + 9 . This would be done 

o 

in the following way: 



p(t) = p(t, 0 Q ) + P ( t , - 9 q ) 



ttY 



^ (rr sinh Y) _1 e 9w (3 a + 8 } 

4*0 o 0 Q -e 

W 

sin J [ir± (0 - 0 J ] 



W 



1-2 exp(~) cos r tt± (9-0 )] + exp (-|^) 

w w "w 



sin J— [*± (6 + 9 ) ] 



“9, 



i 1>U 



-ttY, 



1-2 exp(4_L) cos — [ir±(0 + 0 )1 + exp (-|4i2_) 
0 w w w 



2ttY, 



48 



where each 3 term is actually the sum of two terms, one with 
the term in brackets evaluated using the (+) sign and the 
other using the (-) sign. There appears to be no advantage 
to writing the B-T solution in this way except to show the 
geometrical similarity to Macdonald's result. Indeed, Biot 
and Tolstoy specifically identify both terms as diffraction 
not reflection. True reflections have separately been iden- 
tified by B-T as occurring before the diffraction; an inter- 
pretation that is direct and obvious in their impulse solu- 
tion, but which is obscured in Macdonald's CW solution. 

In summary, the B-T and Macdonald exact solutions have 
been found to be in very close agreement for all of the cases 
examined. Real and image sources have been shown to be a 
valid, if not particularly useful, interpretation of both 
solutions . 
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IV. EXPERIMENTAL APPROACH 



A. GENERAL 

The measurements reported herein were made in the NPS 
anechoic chamber utilizing pulse techniques. The anechoic 
chamber provided the low noise, low reverberation environment 
necessary to make precise amplitude measurements at low source 
levels. The use of a pulsed source provided a signal that 
was both wideband and deterministic and made it possible to 
control the effects of extraneous scattering by judicious 
choice of the sampling window. Detailed descriptions of the 
anechoic chamber and the data acquisition hardware may be 
found in Reference 12, and will not be repeated here, except 
when necessary to explain the results. 

B. DATA ACQUISITION 

A block diagram of the overall experimental data acquisi- 
tion process is shown in Figure 9. Overall timing of the 
acquisition process is controlled by the programmable timing 
simulator. A typical sequence is as follows. The waveform 
generator is triggered by the timing simulator, sending a 
single replica of the preprogrammed pulse to the power ampli- 
fier. This pulse is amplified to approximately 50 volts 
peak, added to the 150 volt D.C. bias voltage, and applied 
to the source (source and receiver characteristics will be 
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Figure 9* Experimental Data Acquisition General Block 
Diagram. 
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discussed separately) . The transmitted signal is received by 
the H inch condenser microphone (B+K 4133) and sent to the 
signal conditioning equipment. Here the signal is amplified, 
filtered to remove low frequency noise and to prevent alias- 
ing in the sampling process, further amplified and digitized 
(12 bits) . Here the timing simulator determines the start 
time and duration of the sampling of each pulse. By selecting 
the sampling frequency and duration, the number of samples 
available to the FFT can be determined. In practice, the 
number of samples is set slightly higher than the number to 
be transformed since the computer begins with the first 
sample and disregards those coming after the specified number 
have been acquired. This process is repeated for each pulse 
with a pause of approximately 80 milliseconds between pulses. 

The computer accumulates the samples by averaging up to 
10,000 samples in the time domain and then computes the 
Fourier transform of this averaged waveform using a standard 
software FFT algorithm. Since the FFT must be taken only 
once, the system runs in real-time, which allows a very 
large number of pulses to be averaged in a reasonable amount 
of time. Since this averaging is coherent, the signal-to- 
noise ratio is proportional to N, the number of pulses, 
rather than the usual N 2 . The averaged time waveform or 
FTT results are stored on a floppy disk for later use. 
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C. SOUND SOURCE 



The ideal sound source for these experiments would be a 
point source capable of high acoustic output over very short 
time durations, with zero output at all other times. These 
characteristics can only be approached with a real source. 
The obvious choice of a spark source is rejected in order to 
avoid finite amplitude effects. Fortunately, the extremely 
low noise levels of the anechoic chamber plus coherent 
averaging makes accurate measurements possible with a source 
of relatively low intensity. The search for a useful source 
of small dimensions and narrow time response has been a con- 
tinuing one at NPS[cf. Ref. 14]. The source chosen for the 
experimental work in this paper were the B § K % inch and 1 
inch condenser microphones. The use of these microphones as 
sound sources is described in Reference 12 and in a B § K 
technical pamphlet [Ref. 15]. In general, the inch source 
covered the higher frequencies (6-40. kHz) while the 1 inch 
source was used to provide additional low frequency coverage 
(down to approximately 1.5 kHz). The 1 inch source is 
actually capable of providing sufficient output over the 
entire frequency range but directivity does become signifi- 
cant at the higher frequencies. 

Typical source waveform and frequency spectra are shown 
in Figures 10 through 13. In Figure 10, the ivaveform at 
the output of the power amplifier is the top curve with the 
received waveform shown below. The received ivaveform (and 
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Figure 10. Typical Time Domain Signal, i" Source. 
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Figure 11. Typical Frequency Spectrum of Acoustic Pulse, i" Source 
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Figure 12. Typical Time Domain Signal, 1" Source. 
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Figure 13. Typical Frequency Spectrum of Received Pulse, 1" Source 



associated frequency spectrum) is actually due to a combina- 
tion of source and receiver response, however, the relative 
contributions of each are not particularly important since 
all results are normalized by an appropriate reference mea- 
surement using the same source and receiver. The time dura- 
tion (or "ringing") of the pulse is important since all of 
the desired diffracted pulse must be sampled before the 
unwanted arrivals from the other edges of plate. The longer 
the acoustic pulse, the larger the specimen must be to 
separate the pulses. This ringing is especially important 
in measuring backscatter since the low amplitude diffracted 
pulse will arrive after the direct/reflected pulse which is 
much higher in amplitude. Figure 11 shows the energy 
spectral density of acoustic pulse in Figure 10. 

Figures 12 and 13 present the same data as Figures 10 
and 11, with the same receiver, but using the 1 inch source. 

In this case, the received pulse was low pass filtered at 
14 kHz in order to enhance the low frequency portion of the 
spectrum. 

Once an acceptable source and receiver were selected, a 
mounting system had to be devised which allowed the micophones 
to be positioned accurately yet present a minimum of 
additional scattering surface. The approach used is shown 
in Figure 14. The diffracting object (in this case a thick 
plate) was suspended with nylon fishing line from a frame 
attached to the ceiling of the achechoic chamber. Heavier 
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Figure 14. Typical Anechoic Chamber Measurement Apparatus. 
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objects were also supported from below. The source and re- 
ceiver were each attached to a pair of thin stainless steel 
wires mounted firmly to the ceiling frame and resiliently 
to the wire mesh floor of the chamber. Source and receiver 
were positioned in the horizontal by moving the ceiling and 
floor attachments and vertically by moving up or down the 
wires. Additional support of the signal cables was provided, 
as necessary, with nylon line from the ceiling frame. 

D. PLATE CONSIDERATIONS 

The plate used in the experiment should be as thin and 
yet, as rigid and non-transparent as possible. In addition, 
it was desirable to keep the overall weight of the plate low 
to facilitate handling and suspension within the anechoic 
chamber. Aluminum plate .475 cm thick was chosen as a 
reasonable compromise. The transmission loss, calculated 
using the approach of Reference 12, was greater than 40 dB 
at 1.5 kHz, increasing with frequency. This \^as considered 
adequate since no thin plate measurements were planned 
where the direct path through the plate was to be substan- 
tially less than the diffracted path. For the thick plate 
tests the direct path \\rould be short but two plates would be 
used and the resulting attenuation would certainly be adequate. 
These estimates were verified by examining the received sig- 
nal prior to the arrival of the first diffraction to ensure 
there was no significant direct path arrival. 
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V. 



BACK SCATTER EXPERIMENTAL RESULTS 



In Figures 5 and 6 it was shoxvn that the B-T and H-K 
solutions give dramatically different results for backscatter 
near the surface of the rigid half-plane. A simple experi- 
ment was performed to determine which solution is correct. 

The one inch B+K microphone was used as the source with the 
one-half inch B+K microphone used as the receiver. They were 
taped together and positioned at a distance of 25 cm from the 
edge at various angles as shown in Figure 15. Although the 
source and receiver acoustic centers are separated by approxi- 
mately 2.5 cm they will be assumed to be coincident. 

The transmitted waveform was the half triangle discussed 
earlier but the received signal was loxv pass filtered at 
60. kHz vice the 14. kHz shown in Figures 12 and 15, to 
achieve a higher frequency response. Each transmitted pulse 
actually resulted in three received pulses. The first pulse, 
which arrived almost immediately, was the direct source/ 
receiver path. The second was the reflection from the plate, 
while the third ivas the desired diffraction from the edge. 
Unfortunately, the first two pulses were considerably higher 
in amplitude than the diffracted pulse. Even with the ability 
to exclude the first two pulses in the time domain, the 
ringing of the reflected pulse interfered with the analysis 
of the diffracted pulse. This limited the geometry to 
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Figure 15. Backscatter Experimental Geometry. 
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source/receiver locations near the plate, which presented no 
particular problem since this was the region of most interest 
anyway . 

The results of this experiment are presented in Figure 16, 
with the free-field pressure at 2r used as the reference. 

For the comparison, the theoretical B-T and H-K results are 
shown for angles approximating those of the measurements. 
(These theoretical results are the same as those presented 
in Figures 5 and 6, now plotted as a function of frequency 
for specific angles.) The B-T theoretical results are 
plotted as a single line since they vary by less than one 
dB over the range of angles presented. I t is clear that the 
measured data tend to agree with the results predicted by the 
B-T solution rather than going to zero at the surface of the 
plate as in the H-K solution. The considerable scatter in 
the measured data, especially near 9 kHz at 30 degrees, is 
attributable to interference between the diffraction pulse 
and the ringing of the reflected pulse. 

The measured data not only favor the B-T solution but 
also confirm the presence of a scattered wave near the sur- 
face. If this scattered wave also exists on the surface, 
as predicted by B-T, (there is no reason to think otherwise) 
it points out a serious fallacy in the Kirchhoff assumption. 

In order to evaluate the integral in Equation 11-25, the 
wave field at all points on the surface of integration had 
to be known. This was accomplished, through the Kirchhoff 
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Figure 16. Comparison of Measured Backscatter from a 
Half-plane with Theoretical B-T and H-K 
Solutions for r = 25 cm at various 9’s (§=© ) 
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assumption, by simply taking the wave field on the surface 
to be the reflection of the incoming wave. This is equiva- 
lent to a boundary condition which forces the diffracted 
wave field due to the edge to zero at the surface of the 
plate. When the solution is transformed to the time domain 
(Equations 1 1 1 - 6 and III-7) where the diffracted and re- 
flected terms are considered separately, the diffracted 
term must go to zero at 9 = 0° in order to satisfy the 
imposed boundary condition. That this is clearly not the 
case indicates that the H-K formulation should not be used 
to estimate the backs cattered wave field near the surface 
of the plate. For forward scattering, the H-K diffraction 
solution would also be expected to go to zero near the Plate. 
Since Bremhorst's data [Ref. 12] show the presence of a 
forward scattered wave at the surface, the above argument 
can probably be applied to this case as well. 
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VI. FINITE PLATE FORWARD SCATTERING EXPERIMENTAL RESULTS 



A. GENERAL 

When a rigorous diffraction solution such as Biot-Tolstoy 
is applied to the problem of diffraction by a real object, 
one must be careful to ensure that the original boundary 
conditions are applicable. For a finite barrier, two impor- 
tant considerations are whether the barrier is large enough 
to be considered infinite in extent and thin enough to ignore 
thickness (the double edge) . Since the geometry of the half- 
plane diffraction problem is such that everything is described 
in terms of a cylindrical coordinate system with the origin 
at the edge, there is a tendency to think of diffraction as 
simply an edge effect. This is the interpretation of the 
H-K boundary diffraction wave. This is also the simplist 
interpretation of a time domain, impulse solution such as 
B-T, where it is easy to visualize an expanding spherical 
wavefront intersecting the edge and reradiating acoustic 
energy as it "propagates" along the edge. If this were 
generally true, one could solve any diffraction problem by 
describing the diffracting object in terms of a series of 
edges (or, more appropriately, wedges) and adding their con- 
tributions. The problem with this approach is that the 
original solution assumed the presence of a rigid surface 

at some angles 0=0 and 0=0 extending to infinity in the 

w 
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radial and axial directions and an infinite fluid elsewhere 
(see Figure 1). When the physical situation approximates this 
condition, the diffraction problem can generally be treated 
by considering a source at the edge. For an extreme example 
where this is not true, one could consider the difference 
between the diffraction due to a rigid half-plane and that 
due to a large plate shrunk to become a thin wire at the 
location of the edge. 

A related prob lem concerns the treatment of barrier cor- 
ners. Here the barrier may be "high" enough to consider the 
radial dimension to be infinite but its width (the axial, or 
z, dimension) is finite. A roadside billboard or a tall 
building could be described this way. If the ground is 
totally absorbant (no reflections) Medwin [Ref. 8] treats 
the diffraction from this type of a barrier by adding the 
impulse response of each of the three edges, truncating the 
solution for each edge when the wave reaches the corner, to 
obtain the total time domain solution. The assumption 
implicit here too is that the boundary conditions of the 
theoretical B-T solution are sufficiently well approximated 
by the actual barrier. 

The third major difference between the theoretical B-T 
solution and the actual barrier is the double edge caused by 
the finite thickness. The theoretical solution is based on 
a rigid, infinitely thin screen whereas the actual barrier 
is often approximately rigid but of finite thickness. 
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Some researchers have minimized this problem in their mea- 
surements by placing a "knife edge" at the top of the barrier 
[ cf . Kawai, Ref. 11] but this introduces additional sources 
of diffraction where the face of the barrier changes angle. 

In any case, the real barrier will always have some finite 
thickness dimension and it is desirable to be able to calcu- 
late or at least understand the nature of its effect on the 
diffracted field. 

A series of experiments was performed in order to define 
these "finite barrier" effects. 

B. FORWARD DIFFRACTION AT A CORNER 

Bremhorst [Ref. 12] showed excellent agreement over a 
wide frequency range between the B-T infinite plate solution 
and pulse measurements on a 5/16" thick plate, provided the 
pulse is sampled in such a way as to exclude late arriving 
diffracted pulses from other edges. However, discrepencies 
between measurements and theory of 5-10 dB were observed at 
frequencies approaching 50 kHz with source and receiver near 
the plate. The purpose of this experiment was to measure 
the diffracted pressure field in the vicinity of the corner 
of a plate and compare it to the field predicted by B-T 
using Medwin's truncation approach. The present geometry 
was chosen to ensure that only the edges intersecting at the 
corner of interest would contribute to the measurements 
(Figure 17). 
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Figure 17. Corner Diffraction Experimental Geometry. 
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The source was placed at the coordinates (x g , y s , z g ) 
of 25.5 cm, 8.0 cm, 8.0 cm, respectively. Separate measure- 
ments were made using the Jg and 1 inch sources to achieve 
the widest possible frequency coverage. The receiver was 
positioned along the x axis (y = 0,z r = 0) at various points 
between r = 3 cm and r = 45 cm. This geometry was dictated 
by two considerations; source/receiver directivity and co- 
herent interference. Source and receiver are both to be 
omnidirectional in the theoretical B-T solution. In order 
to approximate this adequately, it was considered important 
that the source appear omnidirectional when viewed from 
points along a significant portion of the edge near the 
corner. By orienting the source and receiver in this way 
all points on the edge between the corner and the point where 
the least time path crosses the edge were within 2 dB of 
omnidirectionality up to approximately 50 kHz. The symmetry 
of the geometry was chosen to eliminate the frequency domain 
interference pattern created by taking the Fourier transform 
of a time signal which has two coherent pulses separated by 
a time delay. The geometry ensures that the diffracted pulses 
from both edges arrive at the same time. This was accom- 
plished by moving the receiver slightly to maximize the total 
pulse signal at the nearest and farthest ranges and locating 
the intermediate measurement points on a straight line 
between them. 
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These considerations resulted in a geometry which was 
somewhat artificial but would allow a relatively simple analy- 
sis of the corner effect. The Fourier transformed diffrac- 
ted pulse data for the 1" and h” sources are presented in 
Figures 18 and 19 respectively. The free-field pressure at 
a range of 28 cm (the distance from source to corner) was 
used as a reference pressure. Although the measured data 
using both sources show good agreement in the frequency 
region where they overlap, their dimensions are large enough 
to be considered as an effect on the results and are there- 
fore presented separately. The theoretical results from 
truncating and adding the B-T solutions for the two 
perpendicular edges is also shown. The agreement between the 
measurements and this superposition of the truncated 
components of two cases of B-T theory is seen to be excellent 
across the entire frequency region. 

The theoretical and measured results in Figures 18 and 
19 also clearly show the effect of moving closer to the 
shadow boundary. Because of the geometry shown in Figure 17, 
decreasing r also moves the receiver nearer to the shadow 
boundary. The theoretical solution as the receiver approaches 
the shadow, is shown to approach the delta function form 
(Appendix A). The measured data clearly show this result, 
as they approach a flat spectrum for small r, with the more 

_ i, 

familiar f 2 frequency dependence observed at larger r. 
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Figure 18. Forward Diffracted Pressure as a Function of 
Frequency at Various Distances (r) from the 
Corner, 1" Source (Geometry of Fig. 17). 
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Figure 19. Forward Diffracted Pressure as a Function of 
Frequency at Various Distances (r) from the 
Corner, ?" Source (Geometry of Fig. 17). 
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In Figure 20, these same measurements are presented as 
diffracted pressure as a function of range (from the corner) 
for selected frequencies, once again compared to the sum of 
"truncated" B-T theoretical solutions. Also presented for 
comparison are lines of arbitrary level with slopes of -10 
log r (cylindrical) and -20 log r (spherical). An apparent 
source at the corner would be expected to radiate spherically 
while an infinitely long edge would radiate cylindrically . 
Because of the geometry the ranges to the corner and to the 
point on the edge where the least time "ray" crosses are 
essentially equal so the data may be examined directly for 
evidence of a spherical divergence from the corner. 

The data, as well as theoretical results, of Figure 20 
range from a slope of less than 3 dB per double distance 
near the corner, low frequency, to almost 6 dB per double 
distance at the highest frequency and largest distance from 
the corner. In terms of the dimensionless parameter kr 
(wavenumber times range) , it appears that the corner cannot 
be characterised as a simple spherical or cylindrical source 
for low values of kr. As kr is increase, the range depen- 
dence gradually increases to a value approaching 6 dB per 
double distance. Although the results are somewhat incon- 
clusive due to limited data, the two line sources, one at 
each edge, begin to resemble a virtual point source at 
higher values of kr. 
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Figure 20. Forward Diffracted Pressure as a Function of 

Distance from the Corner for Various Frequencies 
(Geometry of Fig. 17). 
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An additional experiment was performed to better define 
the difference between diffraction by an edge and diffraction 
by a corner. The plate used in this experiment was identical 
to that used for the corner (Fig. 17) except that the source 
and receiver were set up on a long edge, as shown in Figure 
21. As before, the measured data are compared to the B-T 
calculated results (for a single infinite edge) in Figures 
22 and 23. Unfortunately, the agreement between measurements 
and theory is not nearly as good in this case as it was for 
the corner. The reason for the significant deviation between 
measurement and theory at higher frequencies and small ranges 
is not known. Figure 24 presents a direct comparison between 
corner and infinite edge results, as a function of range, 
along with the appropriate B-T theoretical results. 

From this comparison it can be seen that the diffracted 
pressure field behind the corner is approximately 6 dB higher 
than that behind a half-plane (single edge) for the larger 
ranges and higher frequencies. As the receiver moves closer 
to the plate, the difference between corner and single-edge 
data decreases. Furthermore, there does appear to be a fre- 
quency dependence in the difference between the corner and 
edge at a given range. To better define this frequency 
dependence and to estimate the effect of moving away from 
the shadow boundary, a final experiment was performed. 

Referring to the geometry of Figure 17, the receiver was 
placed distances (r) of 18 and 32 cm from the corner on a 
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Figure 21. 
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Figure 22. Forward Diffracted Pressure as a Function of 
Frequency at Various Distances (r) from the 
Edge, 1" Source (Geometry of Fig. 21). 
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Figure 23« Forward Diffracted Pressure as a Function of 
Frequency as Various Distances (r) from the 
Edge, i" Source (Geometry of Fig. 21). 
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Figure 2£f. Comparison of Forward Diffracted Pressure as a 
Function of Distance (r) for the Corner and the 
Single Edge for Various Frequencies. 
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line perpendicular to the plate (this is the same geometry as 
the previous corner measurements). Measurements were also made 
at distances of 18 and 32 cm from the corner on a 45° line 
extending into the shadow region "behind" the plate as shown 
in Figure 25. As before, four identical measurements were 
made along a single "infinite" edge for comparison. The 
magnitude of the Fourier transformed data from the corner, 
divided by the corresponding data from the single edge, are 
presented in Figure 26. The upper two curves present data 
taken on the perpendicular line while the lower two are taken 
on the 45° line (deeper in the shadow) . Also shown with the 
upper curves is the theoretical difference between the corner 
and edge based on B-T. The general trend indicated by the 
upper two curves is one of less than a doubling of pressure 
(due to line sources) at the lower frequencies, increasing 
to slightly more than double the pressure at the higher 
frequencies . 

When the receiver is moved deeper into the shadow region, 
the character of the ratio plotted in Figure 26 changes 
dramatically, as shown by the lower two curves. The monotonic 
frequency dependence seen near the shadow boundary has been 
replaced by a more complicated dependence. The average value 
of the difference appears to be approximately 6 dB , even at 
the lower frequencies. The detailed frequency dependence 
will be addressed in the next section. The general conclu- 
sion to be drawn from Figure 26 is that the farther the 
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Figure 25. Experiment Geometry for "Near Shadow Boundary" 
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Near the Shadov; Eoundary (see text). 



83 



receiver moves into the shadow region, the more the diffrac- 
ted pressure field looks like the coherent addition of two 
infinite edges, plus a small frequency dependent variation. 

In summary, the measurements presented herein show that 
the basic approach of adding impulse diffractions from each 
edge in the time domain is valid. Near the corner, the 
diffracted pressure is less than a simple addition of two 
infinite edge diffractions. In this region, the approach 
proposed by Medwin [Ref. 8] , wherein the time domain solu- 
tion is "truncated'' at the point where the edge wave reaches 
the corner, estimates the diffracted pressure very closely 
as shown by Figures 18 and 19. The usual treatment of the 
corner as a source of spherically diverging diffracted sound 
is not evident in the measured or B-T theoretical results, 
at least over the lower range of kr presented. At high 
values of kr , the results may be approaching spherical 
divergence . 

In the high frequency limit, Keller [Ref. 16] predicts 
that the corner does produce spherically diverging waves but 
that the resulting diffracted field decreases with increasing 
frequency much more rapidly than the field produced by an 
edge. This characteristic has been used as justification 
for disregarding the diffraction at a corner when approaching 
practical problems [cf. Ref. 17]. Figure 26 would seem to 
indicate that radial distance from the corner and frequency 
are not sufficient to determine when to disregard the effects 
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of a corner; one must also consider the relationship of the 
receiver location to the shadow boundary. The dimensionless 
parameter implied in Medwin's truncation approach is kd where 
d is the distance from the corner to the point \\rhere the 
least time "ray” intersects the edge. No substantiation 
for this parameter is offered here, however it does take 
into account both distance from the corner and nearness to 
the shadow boundary. 

C. FORWARD DIFFRACTION BY A STRIP 

The next experiment was designed to investigate the be- 
havior of diffraction from an "infinite" strip of material 
whose width was of the order of a wavelength. Two strips 
were actually measured, x^ith \\ridths 10 cm and 4 cm. Given 
the measurement capability of 1.5 to 40 kHz, these two strips 
provided a total range of ka (a represents the width of the 
strip) of 1.09 to 72.8. Source and receiver were located on 
a perpendicular line through the center of each strip, as 
sho\^n in Figure 27a. The source was 14.5 cm from the strip 
for all tests while the receiver was located at distances 
from the strip of 9.5, 23.2, 39.5, and 59.5 cm in the case 
of the 10 cm wide strip and 9.5, 23.2 for the 4 cm strip. 

Once again, this geometry was chosen to eliminate the cohe- 
rent interference that would result from different path 
lengths around each edge. Identical measurements were also 
performed on a single "infinite" edge as shown in Figure 27b. 
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Figure 27. Strip Diffraction Experimental Geometry. 
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As long as linear superposition of edge effects holds, the 
diffracted pressure at any of the measurement points will be 
double for the strip compared to the corresponding plate 
(single edge) , since there is no truncation involved. With 
this in mind, all data are presented as the following ratio: 

p _ diffracted pressure behind strip 
rat diffracted pressure behind plate 

When superposition holds, there will be a simple doubling of 

pressure and this ratio will be +6dB at all frequencies. 

This ratio, calculated from the measurements on the 10 
cm strip and corresponding plate, is presented in Figure 28. 
The data are not presented above 30 kHz because of limita- 
tions in the determination of the perpendicular to the strip 
(i.e. the path lengths around the two edges could not be 
equalized well enough to prevent cancellation effects above 
this frequency) . The oscillating pattern in the low frequen 
cv portion of the measurements is evident from Figure 28. 

If this is interpreted as coherent interference between two 
signals, two features can be extracted from the data. First 
this interference pattern does not depend on distance from 
the strip/plate. Second, by estimating the spacing (in fre- 
quency) between adjacent peaks or between adjacent nulls the 
apparent path difference can be calculated from the relation 
d = \ = 7 Af in meters. From the data of Figure 28, d is 
estimated to be 9 cm, quite close to the strip width of 10 
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(Geometry of Fig. 27). 
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cm. To verify that strip width is the characteristic dimen- 
sion, the experiment x>ras repeated for the 4 cm strip. 

The 4 cm strip results are presented in Figure 29. The 
strip width from Figure 29 data is estimated to be 4.1 cm, 
confirming that the width of the strip is the factor deter- 
mining the form of the "interference pattern". 

The data measured at 9.5 and 23.2 cm are averaged and 
plotted as a function of ka for both strips in Figure 30. 

For simplicity, the data from both 1 inch and 1/2 inch sources 
were averaged in the overlapping region 6-10 kHz. Figure 30 
clearly shows that the frequency dependence is directly re- 
lated to the strip width a, with a ka spacing of 2tt beti^een 
peaks as expected from a simple interference pattern. 

The physical explanation offered for this interference 
is called secondary diffraction or secondary scattering. 

The application of secondary diffraction to the strip is 
illustrated in Figure 51. The figure shows only the 
interaction at one edge; the mirror image of the process 
shown is occurring simultaneously at the other edge. The 
solid lines in Figure 31 depict the path of the pulse as it 
diffracts over the upper edge, while the dashed lines shoiv 
how the diffraction of the same pulse at the opposite edge 
generates signals on both sides of the strip which propagate 
across the strip and diffract once again at the upper edge. 

The doubly diffracted pulse arrives at the receiver at time 
equal to the width ( a ) divided by the propagation velocity 
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Figure 29. Comparison of P , as a Function of Frequency 
at Various Distances (r) for the 4 cm Strip 
(Geometry of Fig. 27). 



90 



10 era 




91 



Figure 30. Comparison of P t as a Function of ka for the 2+ and 10 era Strips. 
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Figure 31* Secondary Diffraction (Dashed Line) Concept Applied to the Strip. 



(c) later than the initial pulse, resulting in the observed 
interference pattern. This late pulse would have two dis- 
tinct characteristics, arising from the additional diffrac- 
tion, which are immediately observable in the measured data 
of Figure 30. First, it will be considerably reduced in 
amplitude relative to the earlier (single diffraction) 
pulse. This is seen in the data as a weak interference 
pattern (the peaks are approximately 6 + 1 dB compared to 
6 + 6 dB if the two were of equal amplitude). Second, the 
difference in amplitude between singly and doubly diffracted 
pulses will be greater with increasing frequency since each 
diffraction should increase the frequency dependence and 

decrease the strength of the diffraction by approximately 
- 1/2 

f . This causes the decaying amplitude of the interference 

pattern with increasing frequency. 

A method of calculating the secondary diffraction has 

been proposed by Medwin^ and implemented by Ms. Emily Childs. 

The general description of the implementation which follows 

2 

is due to Childs. When a spherical pulse intersects an 
edge, a boundary or edge wave is generated which aopears 
to propagate along the edge. Calculating the B-T time domain 
solution in discrete time increments, this edge wave becomes 
a series of discrete sources on the edge, each with a 



^Personal communication to author from Dr. H. Medwin. 
“Personal communication to author from Ms. E. Childs. 
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characteristic amplitude and phase. Extending this to a two 
edge system (the thin edged strip in this case), these dis- 
crete sources on the first edge (now called secondary sources) 
each diffract at the second edge and contribute to the dif- 
fracted field at the receiver. The proper amplitude is 
assigned to each secondary source through a system of vir- 
tual receivers located at the second edge. The phase of 
each source is taken into account by performing the entire 
calculation in the time domain. Figure 52 presents the 

results of this theoretical calculation in terms of P . for 

rat 

each strip at a receiver distance of 23.2 cm., compared to 
the measured data from Figure 30. Two theoretical curves 
are presented for each strip; one which considers the 
secondary diffraction on both sides of the strip (double 
secondary path) and one which considers secondary diffraction 
only on one side (single secondary path). Although the 
double secondary path is the intuitively correct 
interpretation, the measured data obviously agree with the 
single path approach. One possible explanation is that in 
using the double path interpretation, each path should be 
reduced in amplitude by 1/2 (the diffracted pressures along 
the paths on either side of the strip are of equal amplitude 
but opposite sign). Keller [Ref. 16] describes a similar 
phenomenon for the thin slit in an opaque plane using the 
Geometrical Theory of Diffraction. 
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Figure 32. Comparison of Measured and Theoretical P as a Function of ka for 
the if and 10 cm Strips. ra 



At this point it is also appropriate to point out that 
the use of secondary diffraction in the preceding explanation 
of the interference pattern would not be possible with the 
Helmhol tz-Kirchhof f approach, since the scattered or dif- 
fracted field on the plate would be zero. There would be no 
way the pulse could propagate between edges. 

It appears that the secondary diffraction explanation can 
also be applied to the somewhat more complicated geometry of 
the corner. Recalling Figure 26, there was a detailed fre- 
quency dependence in the results, especially where the 
receiver was deep in the shadow region (the lower curves of 
Figure 26). Consider now a corner with one edge horizontal 
and one vertical as shown in Figure 17. There will be a pulse 
that travels from source to receiver by diffracting over the 
horizontal edge. In addition, a pulse also diffracts at the 
vertical edge, propagates across the plate near the corner, 
and diffracts over the horizontal edge to the receiver. 

Unlike the strip, the difference in path length of the two 
pulses is a function of distance from the corner, an effect 
which can be observed as a difference in the interference 
patterns in the two lower curves of Figure 26. 

D. FORWARD DIFFRACTIONS BY A THICK BARRIER 

The usual theoretical approach to solving the problem 
of diffraction by a plate is to consider it to be an 
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infinitely thin, rigid half-plane. However, a real barrier 
which significantly attenuates the direct path of sound must, 
of necessity, have a thickness which is a significant frac- 
tion of a wavelength at high frequency. In Reference 12, 
Bremhorst presented data which showed the diffracted pressure 
at various points in the shadow of a 3/16” thick steel plate. 
In particular, Figures 37 through 46 of Reference 12 showed 
that the measured diffracted pressure was generally less 
than that predicted by B-T, especially near the surface of 
the plate (cm the shadow side) . Bremhorst attributed this 
to the finite thickness of the plate. To confirm this 
hypothesis, an experiment was performed using two barriers 
of different thickness. 

The variable thickness was achieved by placing a piece 
of aluminum stock of the appropriate thickness between two 
3/16” aluminum plates and machining the resulting "edge” 
smooth and square. Damping material was placed in the gao 
between the plates to attenuate the direct path. An example 
of the thick edge can be seen in the photograph of Figure 
14. Thicknesses of 2.2 and 3.5 cm were evaluated utilizing 
the geometry shown in Figure 33. Measurements were made at 
9 = 300° and 345°, where there was considerable disagreement 
between Bremhorst's measurements and B-T theory. The measured 
data are shown in Figure 34, compared to the Bremhorst's 
results and the B-T theory. 
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Figure 




< — t — > 



33 * Thick Barrier Measurement Geometry. 



98 



Pressure in dB re free- field at 50 cm 



Measured Data: 



^ t- .48 cm (Bremhorst) 





I I I I 

5.0 10.0 20.0 50.0 

Frequency in kHz 

Figure 34. Forward Diffracted Pressure as a Function of 
Frequency for Various Thickness (t) Plates, 
© 0 = 15 °, r s r = 25 cm. 
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The data presented in Figure 54 indicate that barrier 



thickness is an important parameter in determining the dif- 
fracted pressure field, especially when the thickness is on 
the order of one wavelength (10-15 kHz for the thick barriers 
tested) . It has been suggested by Medwin ^ that this thick- 
ness can be handled by a secondary diffraction calculation 
based on the B-T theory, as discussed earlier for the thin 
strip, by considering the thickness to be a double edge. 

This has not been verified at this time, although the thick- 
ness does appear to explain the discrepancy between Bremhorst's 
data and the B-T theorv. 



E. FORWARD DIFFRACTION BY A BARRIER ON A RIGID BASE 

A typical noise control barrier consists of a plane 
screen or plate mounted perpendicular to a flat base, as 
shown in Figure 35. The total pressure field at the receiver, 
if it is in the shadow region, consists of the source signal 
arriving via multiple paths. The usual approach to solving 
this problem is to consider that for the top edge there are 
four geometrically distinct paths, each involving one 
diffraction at the top of the barrier, made up of all possible 
combinations of source, receiver, and images of source and 



^Personal communication to the author by Dr . H. Medwin. 
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receiver in the bas e , as shown in Figure 35. Using the no- 



tation of Figure 35, these paths are denoted as SDR, S'DR, 
SDR', S'DR'. It should be pointed out that S' and R' here 
are true images as used in optics and acoustics to describe 
reflection from a plane surface. The term "image" is also 
used by authors in the noise control field [cf. Ref. 5] in 
conjunction with the "image of the source in the barrier", 
as discussed in section II. This usage is considered here 
to be inappropriate and will always be set apart by quotation 
marks ("image") to prevent confusion with the more common 
term (image) . 

To test the hypothesis that the four paths described 
above adequately characterize propagation from source to 
receiver in this case, a simple experiment was performed. 

A barrier with base was constructed of 3/16 inch aluminum 
plate with an overall x^idth of 120cm and barrier height of 
25 cm in the basic configuration shown in Figure 35. Source 
and receiver (both 1/2" B+K microphones) were positioned so 
that each of the four paths would have a different length, 
and each path would intersect the axis of the microphone 
at approximately the same angle. This would alloiv at least 
a qualitative comparison of the relative strengths of the 
various paths. Use of the pulsed source of earlier experi- 
ments would allo\^ the identification, in the time domain, 
of each path by its characteristic propagation time. 
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Similarly , undesirable paths such, as backscatter from the 
edges of the finite base and diffraction around barrier sides 
would have propagation times longer than the longest (S'DR') 
path of interest and could be eliminated in the time domain. 
The results of this experiment are shown in Figure 36, 
compared to the calculated arrival times of each of the four 
paths. It is clear from these results that the four paths 
constructed from the source, receiver and their images in 
the base are the dominant paths in this basic barrier confi- 
guration. If any other paths exist, they must be of rela- 
tively low level or coincide with one or more of the four 
basic paths. The path due to the "image of the source in 
the barrier" as discussed by Isei [Ref. 5] would coincide 
with the SDR and SDR' paths of Figure 36. The relative 
amplitude of the various paths also behave as expected, with 
the later arrivals (longer paths) having lower amplitudes. 

The strengths of the various paths were examined by 
varying the geometry slightly. By increasing the barrier 
height to 42 cm and changing the receiver location (Figure 
37) , the separation in time between the paths was increased 
engough to allow a separate spectrum analysis of each pulse. 
The data of Figure 38, taken with the source/receiver oriented 
to insonify/receive both SDR and S'DR paths equally 
(approximately), show the separation between the two pulses, 
and the general trend in their amplitudes. As expected, the 
S'DR pulse appears to be an exact replica of the SDR pulse. 
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Amplitude 



S' DR' 



SDR S' DR SDR' 




Figure 36. Total Received Time Domain Signal Due to the 
Four Paths of Barrier on Base (Geometry of 
Fig. 35) . 
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Figure 37. Geometry of Barrier on Rigid Base, Tall Barrier 



Amplitude 



NOTE: Relative Amplitudes 




Figure 33. 



Received Time Domain Signal Due to the SDR 
and S’ DR paths (Geometry of Fig. 37). 
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The actual strength of each path was. measured hy "pointing" 
the source in the direction of the path of interest (.the re- 
ceiver always pointed toward the edge, D) and Fourier trans- 
forming only the received pulse corresponding to that path. 
This was done for both the SDR and S' DR paths. The data, 
normalized by the free field pressure at the receiver with 
the barrier and base removed, (source/receiver pointed at 
each other) are presented in Figure 39. Also presented for 
comparison are the B-T infinite half-plane theoretical results 
for each path. The agreement between theory and measurements 
shown in Figure 39 is excellent, indicating that the diffrac- 
ted field at the receiver due to SDR and S'DR paths may be 
computed separately, using the geometry of the source and the 
image of the source in the base. The total diffracted 
pressure at the receiver may then be obtained by adding the 
contributions of each path. 

It is important to note that the theoretical results in 
Figure 39 were obtained from the complete form of the B-T 
infinite half-plane solution. The complete theoretical 
solution for the diffracted field at R in Figure 35 consists 
of the sum of four terms. Isei [Ref. 5] approaches the same 
problem by adding the contributions of six different paths, 
as shown in Figure 3 of Reference 5. It appears that he 
has used the solution due to Macdonald in a form which is 
separated into "real" and "image" components as discussed in 
section II and is applying these separately to the six paths. 
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Pressure in dB re Direct (free- field at 52.3 cm) 
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Figure 39. 



Measured data 
B-T Theory 



I I I J 

5. 10. 20. 50. 

Frequency in kHz 

Diffracted Pressure as a Function of Frequency 
for the SDR and S’ DR Paths Comnared to the 
B-T Theory (Geometry of Fig. 37). 
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Figure 39 shows that the correct approach is to apply the 
complete solution (e.g. in Macdonald’s solution this would 
be the sum of the "real" and "image" source components) to 
each of the four paths to obtain the correct diffracted field. 
One need not consider the "image of the source in the barrier" 
because it is an integral part of the complete solution. 

In addition to the four paths described earlier, it is 
possible to have an additional path which propagates from 
the source, diffracts at the intersection of the barrier and 
base B, then diffracts over the top edge of the barrier to 
the receiver, shown in Figure 40 as SBDR. This assumes a 
non-perpendicular intersection of barrier and base since 
Tolstoy [Ref. 7] shows that there will be no_ diffraction from 
a 90° "interior" corner. The theoretical prediction of this 
diffracted pressure at the receiver was approached in the 
same way as for the thin strip, discussed earlier. 

The geometry shown in Figure 40 was constructed of 3/16" 
aluminum, as before, with the source and receiver located as 
shown. The SDR and SBDR paths [see Figure 40) can be seen 
directly in the time domain in Figure 41. (The relative 
amplitudes are only approximate due to source directivity.) 
Once again, the source was pointed toward the path of interest 
and the received pulse was analyzed separately. The frequency 
spectrum of the SBDR pulse, normalized by the direct, free- 
field pressure, is shown in Figure 42. Also shown is the 
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Figure Zfl. Received Time Domain Signal Due to the SDR 
and SBDR Paths (Geometry of Fig. 40). 
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Figure 42. Diffracted Pressure as a Function of Frequency 
for the SDR and SBDR Paths Compared to the B-T 
Theory and Secondary Diffraction Theory, 
respectively (Geometry of Fig. 40). 
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theoretical result from the secondary diffraction program. 

For comparison, the SDR path data and theory from Figure 59 
are also presented. (The SDR geometry is identical to the 
previous experiment.) The additional diffraction in the 
SBDR path is seen in the data as a different frequency depen- 
dence as well as a lower overall level. When the B-T theory 
is used as a basis for secondary diffraction, there is 
excellent agreement with the measured data. 
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VII. CONCLUSIONS 



The exact impulse diffraction solution of Biot and 
Tolstoy C B - T ) has been compared to the exact CW diffrac- 
tion solution of Macdonald for forward scattering and to 
the impulse solution to the Helmholtz equation (H-K) for 
backscatter using the Kirchhoff approximation. For back- 
scatter from an infinite edge, B-T and H-K have been found 
to reduce to a similar although analytically different time 
domain form. Comparison of B-T and H-K in the frequency 
domain showed substantial disagreement, except for geometries 
where the diffracted pulse was close to the reflected pulse 
in time, where both solutions approached the same value. 
Backscatter near the plate was found to be grossly under- 
estimated by the H-K solution, due to the boundary condition 
used by Kirchhoff to simplify the integral. The backscattered 
pressure near the surface predicted by B-T was confirmed by 
laboratory measurement. 

Macdonald and B-T solutions for forward diffraction from 
a plate were both found to be in good agreement with the 
experiment. Both solutions can be expressed as the sum of 
a term containing the source/receiver angles as ( 0 - 0 Q ) and a 
term containing (9 + 9 0 ). These terms were originally identi- 
fied by Macdonald as the "real" and "image" source terms. 

The concept of "image of the source in the barrier", as used 
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in the noise control literature, is found to be a mis- 
interpretation of Macdonald's solution. For a barrier on 
a rigid base, the diffracted pressure due to the real source 
and the true image source (the image of the source in the 
base) were measured separately and both were found to agree 
with the complete (both 8‘® 0 and 0 + 8 o terms) B-T solution. 

The solution of finite barrier problems was found to 
be facilitated by the time domain, impulse nature of the B-T 
solution. Measurements of forward diffraction near a plate 
corner show that Medwin's proposed approach of truncating 
the B-T time domain solution to account for the corner accu- 
rately predicts the diffracted field. Experimental results 
do not agree with simple spherical radiation from the corner 
over most of the range of kr between 1 and 150. Limited 
data do show a trend toward spherical divergence as kr 
approaches ISO. Measurements of the forward diffraction by 
a thin strip confirm the presence of secondary diffraction, 
i.e., waves that first diffract at one edge, propagate across 
the surface of the strip, and diffract again at the other 
edge. The time domain B-T solution is shown to be a useful 
starting point in characterizing this secondary diffraction. 
At high frequencies, forward diffraction loss by a plate of 
finite thickness is found to be considerably greater than 
the theoretical prediction, for an infinitely thin screen. 

It appears that this may also be explained by the double 
edge creating secondary diffraction with an additional 
diffraction contribution at the second edge. 
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APPENDIX: A 



APPROXIMATION OF THE BIOT -TOLSTOY SOLUTION 
NEAR THE SHADOW BOUNDARY 

Since the B-T solution is a delta function impulse solu- 
tion, it is expected that it will approach the delta function 
form for limiting cases where the receiver is approaching 
the boundary between "illuminated" and "shadow" regions, 
herein referred to as the shadow boundary. As this limit is 
approached, the early time portion of the diffracted pulse 
will increasingly dominate the solution, creating a difficult 
numerical calculation problem. To better understand this 
problem, the solution will be approximated for the following 
conditions : 

A. 9 = 2tt (thin plate) 

W v * ' 



B. 9 = 9 



o 



+ it + e (receiver in the shadow region, e 
small) 



C. t - x -*-0 



o 



from the positive side (early 
diffracted pulse) 



For 9 



w 



2m, the complete solution can be written as 



P(t) 




(A - 1) 



where 




(A - 2) 



and 3 is a four term sum written symbolically as 
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8 



1 



in [ 2 U±0±0o)] 



- 2exp ( -jY) cos [y(7T±0±0 o ) ] + exp(-Y) 



(A- 3) 



For the early diffracted pulse 

At = t - x -*-0 
o 

Medwin [Ref. 8] shows that when ^-<<1, the factor containing 

o 

sinh Y can be approximated as 

-i -I 

(rr sinh Y) = (2 t c 2 rr ) 2 
o oo 

(This same result for the specific case of backscatter was 
shown in section III.) This factor goes to infinity as At 
approaches zero and is generally assumed to be the factor 
that determines the time dependence of p(t) under these 
circumstances. Time dependence also enters A-l through the 
exp(-Y/2) term but this will be ignored because this function 
only varies between 1 and zero and is relatively constant 
for Y (and therefore At) approaching zero. 

The 8 term contains both time and geometry dependence in 
a rather complex way. The approach here will be to separate 
the 8 into four terms and examine each one in the limit as 
both Y and z approach zero. (y is directly, although not 
linearly, related to At. The expansion of arg cosh is not 
trivial and so the limiting procedure here will be with 
respect to Y rather than t.) The 8 term can be written as 
(8 ++ + 8 + _ + S + + 8 _) with the subscripts referring to the 
various combinations as shown in Table (A-l). 
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TABLE A- 1 



TERM 



COMBINATION OF 
TT , 0 , 9 

0 



SUBSTITUTION OF 

9 — 9 + TT + £ 

0 



+ + 



h- 



- + 



TT + 9 + 0 



TT + 9 - 9 



TT - 9 + 9 



2 tt + 29 + e 

o 

2tt + e 

-£ 



B 



TT - 9 - 9 

0 




e) 
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Using the usual, first order approximations for sine and 
cosine of small angles, these terms can be approximated by 
the following 

sin 9 + e/2 cos 9 

a ~ . o o 

^ + + “ 

1+2 ext>(-^/2) [cos9 q - e /2 sin 6 q ] + exp(-Y) 

£ / 2 

8 , ? 

1 + 2exp(-^"/2) + exp(-Y) (A-4) 

e/ 7 

B » 

1 - 2 exp ( - Y/ 2) + exp(-Y) 

g ~ . sin 9 0 + z _j 2 cos 9 0 

1 - 2exp ( -Y/ 2) (cos 9 q + e /2sin 9 q ) + exp(-Y) 

With these first order approximations, the behavior of 
each 8 term can be examined at Y = 0 



„ _ sin Qp + £ /2 cos 9 n 

Y=0 2(1 + cos 9 - 2 sin 9 ) 

„ _ sin + Z J_ 2 cos 9 q 

Y=0 2(1 - cos 9 - z /2 sin 9 ) 

v o' o 





£ 

s 

(denominator goes to zero) 



From this first order analysis, it is seen that three 
of the four 8 terms have the potential of dominating the 
t>(t) solution at small values of Y. These terms will now be 
examined to second order in the denominator. (The numerator 
contains no difference terms and so the first order 
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approximation for small e will be retained for simplicity.) 
Considering 3_ + , the following approximations are used 
cos C ~ £ / 2) = 1 - 

Y 2 

exp (-Y) r 1 - Y + — 



2e 






to obtain 
3 , * 



Neglecting terms higher than second order in the small 
quantities Y and e, 

n „ 2e 



- + 



Y 2 + £ 2 (l-Y/2) 



(A - 5) 



which approaches a finite value as Y approaches 0, 



lim 3 
Y-*0 



- + 



2 

£ 



Using an identical approach yields the following results 



for B ++ and 3 



lim 3 
Y+0 



sin 9 + Z 1 2 cos 9 

o o 



(A- 6) 



+ + 



2(l+cos 9 + /8 cos 9 - £ /2sin 9 ) 

^ o o o 



sm 9 + £/2 cos 9 

1 lm 3 _ = - 

Y+0 ” 2 ( 1- cos 9 o + e 2 / 8cos 9 0 ' e/2sin 9 q ) (A-7) 

A detailed analysis of the 3 ++ and 3__ terms was not 
done, however, they assume the following maximum values at 
Y = 0: 
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3 + + = f at 0 o -• TT ; Y = 0 
B__ = - \ at 9 q = 0; Y = 0 

Since these limiting values are the same as for the 3 term 

and will be less for other values of 0 , it is concluded that 

o 

the S_ + term will dominate 8 for the early diffracted pulse. 

_ - 1 

For comparison, the (rr^sinh Y) term can also be 
approximated for Y approaching zero as 

- 1 Y 3 

(rr sinh Y) = — ^-(Y + r + ) _1 s — - — (A-8) 

o J rr'- 6... v 

o rr^Y 

o 

Comparing Equations A-8 and A.- 5 leads to the conclusion that 
the rr^sinh Y factor will determine the time dependence of 
the early diffracted pulse only when e 2 is much greater than 
Y 2 and Y is much less than 1. 

In terms of the basic diffraction problem geometry, this 
means that as the receiver approaches the shadow boundary 
from inside the shadow (e approaches zero) , the time region 

. i 

over which the (rr sinh Y) factor will determine the time 

o 

dependence of p(t) becomes very small. This region must be 
determined, based on the geometry, before an analytical plus 
digital approach such as outlined by Medwin [Ref. 8] can be 
employed. This is illustrated in Figures A-l, A- 2 whicn show the 

behavior of (sinh Y) 1 and 8_ + as a function of Y for two 
values of e. Since the B-T solution contains the product of 
these terms, one can be assumed to be the factor which deter- 
mines the time dependence only when the slope of the other 
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s' 

s' 




Figure .1-1 Behavior of the B-T Solution Near the Shodow 
Boundary, £-18 . 
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s 




Figure A- 2 Behavior of the Q B-T Solution Near the Shadow 
Boundary, e - 



123 



approaches zero. Based on this, the (sinh Y) 
dominate below approximately Y = .1 for e = 1 
Y = ,01 for e = 4.2°. 



term would 
8° and below 
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